Density wave instability in a 2D dipolar Fermi gas 
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We consider a uniform dipolar Fermi gas in two-dimensions (2D) where the dipole moments of 
fermions are aligned by an orientable external field. We obtain the ground state of the gas in Hartree- 
Fock approximation and investigate RPA stability against density fluctuations of finite momentum. 
It is shown that the density wave instability takes place in a broad region where the system is stable 
against collapse. We also find that the critical temperature can be a significant fraction of Fermi 
temperature for a realistic system of polar molecules. 
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I. INTRODUCTION 

The realization of ultracold molecules confers remark- 
able opportunities to study new states of quantum mat- 
ter, and is also of much interest for quantum comput- 
ing [ij. In particular, the investigation of polar molecules 
is a good candidate for studying various quantum many- 
body states because the anisotropic and long-range na- 
ture of the dipole-dipole interaction offers r|ch properties 
that do not occur in non-dipolar systems [2j . 

Recently K.-K. Ni et al. succeeded in the creation of a 
dense gas of *°K-®''Rb polar molecules with temperature 
T fa 2.5Tp where Tp is the Fermi temperature. For 
the absolute rovibrational ground state of the fermionic 
molecule attained in the experiment, the magnitude of 
the electric dipole moment is 0.566 Debye, yielding to 
significant interaction effects in a degenerate gas of dipo- 
lar fermions at lower temperatures. As shown in Ref. [J], 
the anisotropy of the dipole-dipole interaction results in 
a deformed Fermi surface of the ground state of a dipo- 
lar Fermi gas. Since the Fermi surface is the determining 
factor in low energy properties of the Fermi system, theo- 
retical studies of collective oscillations dil, stability 0, 
expansion dynamics after turning off the trapping po- 
tential [a, B, zero sound propagation in a homogeneous 
system l9l. llO|, a nd equilibrium properties at finite tem- 
peratures [llMlSf under the deformed Fermi surface have 
been done recently. 

Some of most fascinating challenges of a dipolar Fermi 
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gas concern phase transitions of possible ordered phases. 
So far, the realization of superfluid phase in 3D [13, Ull 
and 2D biaxial nematic phase [l^, and ferromag- 

netic phase of two component mixtures [l9| have been 
proposed theoretically. Moreover, there will be another 
interesting ordered phase arisen from the repulsive part 
of the long-range dipole-dipole force, that is a density 
wave phase. 

In a density wave phase, translational invariance is bro- 
ken spontaneously, resulting in a distorted density dis- 
tribution [23] ■ The phase transition of the charge den- 
sity wave in ID conductor was originally discussed by 
Peierls [2l|, where he found the instability of a metallic 
state at zero temperature. In mid 1970s, the charge and 
spin density waves were discovered in experiments and 
a number of interesting phenomena were found in static 
and dynamic properties of the density waves p^ j. 

In the present paper, we consider a uniform dipolar 
Fermi gas in 2D with the dipole moments aligned to 
an orientable external field. We study a realization of 
the phase transition of a density wave by evaluating the 
stability condition for the thermal equilibrium state in 
random phase approximation (RPA). By calculating the 
critical temperature Tc, we will show that the density 
wave can be achieved at a significant fraction of Tp for 
a realistic system of polar molecules. Here we note that 
the density wave is brought purely by the interaction ef- 
fect in contrast to the density waves in an optical lattice 
potential as discussed in Refs. [23l. [24|. 

This paper is organized as follows. Section II presents 
the mean field model of the uniform gas of fermions inter- 
acting via dipole-dipole forces in 2D. Section III presents 
the ground state properties at zero temperature and dis- 
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cusses RPA stability. We obtain a phase diagram at zero 
temperature as a variance of the magnitude of the in- 
teraction strength and the angle between the direction 
of the dipole moments and normal direction of the 2D- 
plane. Section IV applies our model into dipolar Fermi 
gases at finite temperatures, which leads to the critical 
temperature of the density wave phase transition. Sec- 
tion V is a summary. 



II. DIPOLAR FERMI GAS IN 2D 

We consider a gas of dipolar fermions of mass m 
and electric or magnetic dipole moment d. The dipoles 
are confined by a harmonic trapping potential V{z) — 
with a trap frequency Wz in the z-direction. 
For hu)z ^ ep where ej? is the Fermi energy, the system 
is effectively 2D. The dipole moments are alined by an 
external electric or magnetic field, E, subtending an an- 
gle 9o with respect to the z-axis as shown in Fig. [1] In 
this case, the Hamiltonian of the system is given by 

i=i t^j 

where N is the number of fermions and is the position 
vector of «th-particle in the x — y plane. The 2D dipole- 
dipole interaction in Eq. ([T]) is described by 



Vddir) 
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{l — 3 sin^ ^0 cos^ ip} 
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= — jP2(cos0o)- 2^"^ (2) 



where is the azimuthal angle relative to the a;-axis, see 
Fig.[TJ In the above expression we have 2nd-order Legen- 
dre polynomial, P2{cos9o) = (3cos^ 9o — l)/2. According 
to the procedure introduced in Ref. [lOj, the momen- 
tum representation of the 2D dipole-dipole interaction of 
Eq. ^ is given by Vdd{q) ^ Vq + V2d{q) where 

Vq = 27rd2p2(cos0o)-, (3) 

rc 

V2diq) = 7rd\{-2P2(cos6io) +sin2 6locos20} . (4) 

In Eq. (|4]), is the angle with respect to the g^^-axis. In 
Eq. ([3]), rc is a cut-off length and an order of size of a 
dipolar particle, say size of a polar molecule, where the 
dipole moment is no longer ideal dipole moment. Thus 
the divergence in the limit of Tc — ?► is artificial [1^ . As 
we shall see below, however, the total energy for dipolar 
fermions is always finite even in the limit of rc — > 0. 
Then we do not need to pay attention to such an artificial 
divergence. 

The thermal equilibrium state of dipolar fermions at a 
temperature T is assumed to be translational invariant 
in the x-y plane. In Hartree-Fock (HE) approximation. 
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FIG. 1: (Color online) Experimental setup of 2D dipolar 
fermions. The dipole moments are alined by an orientable 
external field, E, that forms the angle Oq with respect to the 
z-axis. ip denotes the azimuthal angle with respective to the 
a;-axis of the relative coordinate of two particles in the x-y 
plane. 



the total energy per area derived from the Hamiltonian 
of Eq. ([!]) is represented by 



Etot — Ekin + Eint , 
Ekin — 
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■f{k)f{k')V2d{k-k%7) 



2 J (27r)2 J (27r)2 
where f{k) is the Fermi distribution function at T 



f{k) 



1 



(8) 



exp{(e(fc)-^)/fcBr} + l' 

with Boltzman constant fc^ and the chemical potential 
/i. The single particle energy e(fc) is obtained by solving 
the selfconsistent HE equation 

<k) I ^V2d{k - k')f{k'), (9) 



2m 



(27r) 



under the constraint for the number density of fermions 
in 2D 

We note that the mean field interaction energy Eint 
of Eq. ([7]) is the sum of the direct and exchange energy, 
where a total cancellation of the short range interaction 
Vq of Eq. ([3]) occurs due to Fermi statistics. Thus the 
divergence of Vq in the limit of rc — >■ disappears in E^t 
or Etot- 



III. GROUND STATE AND STABILITY AT 
ZERO TEMPERATURE 

In this section, we obtain the ground state at zero tem- 
perature T = in HE approximation where the distribu- 
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tion function is replaced by 



f{k)^Q{eF^e{k)) 



(11) 



with ^ = cp- In the above equation 9() denotes Heav- 
iside's step function. When the external field tilts with 
respect to the z-axis, the anisotropy of V2d{Q) in the 2D 
momentum space causes the deformed Fermi surface. Be- 
fore the study of numerical calculations in HF approxima- 
tion, we introduce a variational approach that captures 
physical insight of the ground state. 



A. variational method 

In Refs. [1,[1]> we developed a variational method (VM) 
to a 3D dipolar Fermi gas, which describes the deformed 
Fermi surface with an ellipsoidal shape. In Ref. [l6| . 
Bruun et al. applied VM to a 2D dipolar Fermi gas with 
the variational density distribution 




-1.5 -1 -0.5 0.5 1 1.5 

kx/kpo 



FIG. 2: Deformed Fermi surface for g = 1.0 and 9o = 
arccos(l/-\/3) derived from numerical calculations in HF ap- 
proximation (cross) and VM (solid line). 
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B. Ground state properties 



where the positive parameter a describes an elliptical 
Fermi surface and kpo = ^/Awri^. When 6*0 ^ 0, the mi- 
nus sign of the interaction energy ([7]) and the anisotropy 
of V2d{(l) tend to stretch the Fermi surface along the x- 
axis leading to a < 1. 

Under the assumption of the distribution function of 
Eq. (fT2|) . the total energy in VM in units of h?n\^/m is 
given by, see Appendix A, 

ftot = |(^ + a')-§5/(a;^o), (13) 

where g = AmcP kpQ / ^ttH^ and I{a; Bq) is defined by 

I{a; Oq) = E{1 - a ) 

a 

+^ |-T^^(i - + ^^(1 - "^) j(i4) 

a [ 1 — 1 — J 

for a < 1 and I{a;9a) = — tt for a = 1. In Eq. ([H]). 
K{m) and E{m) arc the complete elliptic integrals of the 
first and second kind defined by 

K{m) = / dO (15) 

.7r/2 

E{m) = / rfeVl-msin^ei, (16) 
Ja 

respectively. We find the VM ground state by minimizing 
the total energy (|13|) giving rise to an optimized varia- 
tional parameter a = a^. 



Figure m shows the deformed Fermi surface in the 2D 
momentum space in units of kpo for g = 1.0 and 9q = 
arccos(l/\/3) where P2(cos0o) = 0. The cross and solid 
line correspond to results from numerical calculations in 
HF approximation and VM, respectively. It is shown that 
two results match very well as the case for a 3D dipolar 
Fermi gas [9j and the Fermi wave number derived from 
tp = ^{k), kp{(l)), is dependent on 0. In Fig. [Sj we plot 
the aspect ratio of the Fermi surface , kp{n/2)/kp{Q), 
of the HF (cross) and VM (solid line) ground state as a 
function of g for — O.Itt, 0.27r, arccos(l/-\/3). We note 
that the aspect ratio in VM is given by kp(TT/2)/kp{0) = 
a^, see Eq. p2)) . This result reveals that the shape of 
Fermi surface is well approximated by the elliptical shape 
for different values of g and 6*0. 

We numerically calculate the Fermi energy, ep, by a 
partial derivative of the total energy of the system with 
respect to the number density n2d- Figure [4] shows a plot 
of ep in units of h^n2d/'m of the HF (cross) and VM (solid 
line) ground state as a function of g for different values of 
^0- Again two of results agree very well. As g increases, 
ep increases for Oq = O.Itt, 0.27r and decreases for 6*0 = 
arccos(l/'\/3). This result reveals that an average effect 
of V2d{q) turns to be from repulsive to attractive as the 
direction of the external field is further apart from the 
z-axis. 



C. Stability condition 

We examine the stability of the ground state and ob- 
tain a phase diagram at zero temperature in the g-0o 
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FIG. 3: Aspect ratio of the Fermi surface, kF{n/2)/kF{0), as a 
function of g for 6*0 = O.Itt, 6o — 0.2tv, and do = arccos(l/\/3) 
from top to bottom. The solid line, dashed line, and dotted 
line are results of the HF ground state, and the cross, blank 
square, and filled square are results of the VM ground state. 



FIG. 5: Single particle energies e{k, = 0) and e(fc, — 7r/2) 
for g — 1.0 and for 6^0 = arccos(l/v'3) derived from numerical 
calculations in HF approximation (crosses) and VM (solid and 
dashed lines). Dotted line represents the Fermi energy in HF 
approximation. 
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FIG. 4: Fermi energy as a function of g for 9o — O.Itt, 9o = 
0.27r, and do ~ arccos(l/\/3) from top to bottom. The solid 
line, dashed line, and dotted line are results of the HF ground 
state, and the cross, blank square, and filled square are results 
of the VM ground state. 



plane. O wing 
take RPA [m 



to judge the stability of the system, we 
[27I and obtain the stability condition 



1 + V24q)x{q) > o, 



(17) 



where x{q) is the density-density response function 
against density fluctuations of finite momentum q de- 
fined by 



xiQ) 



d^k f{k + q)-f{k) 



(27r)2 e(fe)-e(fc + q) 



(18) 



As we have seen in the previous subsection, the as- 
pect ratio and the Fermi energy in VM approximate to 
those in HF approximation quite well. Thus we calcu- 
late the response function by use of the single particle 
energy that reproduces the variational distribution func- 



tion The single particle energy in VM reads 



e(fc) 
e(0) 
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(19) 
(20) 



where A represents the curvature of the single particle 
energy and is determined by the relationship 



ep = e(0) + A' 



2 1,2 



h k 



2 m 



(21) 



Figure [5] shows the single particle energy in HF approx- 
imation and VM as a function of fc = |fc| at = and 
tt/2. The single particle energy in VM matches with that 
in HF approximation below and near the Fermi energy. 

By use of the density distribution and single par- 
ticle energy p9)) in Eq. (flSl) . the analytic form of the 
response function for q — (q, 0) is obtained by 



x(g) 



27rA2;i 



2fi2 



1- Wi- 



2fcj^(0) 



Q{q~2kFm 



(22) 

where the angle dependent Fermi wave number is given 
by 



k 



FO 



2 r'nG2 /A _1_ _1_ slu'^ 



(23) 



Q!n COS 



We plot x{q) of Eq. ([22)) in Fig. |6] that is constant below 
q < 2kF{4>) and monotonically decreasing for q > 2kF{(j)), 
showing a singular behavior at g = 2kp{(f)) 

Apart from 9q — 0, the 2D dipolar interaction V2d{Q) 
is anisotropic and is most negative at (j) = tt/2 for a 
fixed magnitude of momentum q. Thus the system is 
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FIG. 6: Density-density response function in VM. 



FIG. 7: Phase diagram at zero temperature in terms of g and 



expected to be unstable against density fluctuations of 
the momentum q = (g, 7r/2). Once the condition 



< 



(24) 



is fulfilled, the system becomes unstable. The mini- 
mum value of the left hand side of Eq. occurs at 
q = 2A:F(7r/2) — 2aokpo Thus we obtain the instability 
condition of the system as 



gcos^^o^ > — -0.212. 



(25) 



K "'^ < where 



7 dStot 



1 jt^^^tot 

8^ 



(26) 



This result reveals that the density wave phase transition 
takes place in a broad region for larger g and smaller Oq 
where the system is stable against collapse. 



IV. CRITICAL TEMPERATURE OF THE 
DENSITY WAVE PHASE TRANSITION 



In case the above condition is satisfied, the density fiuc- 
tuation of the momentum q = 2aokp along the y-axis 
starts to develop. This indicates that the phase transi- 
tion from a normal phase gas into a density wave phase 
takes place for larger g and smaller Oq. Note that the 
expected density wave is a planar wave in the transverse 
direction to the dipolar direction. As shown in the insta- 
bility condition (pS)) . the system becomes unstable even 
when ^0 = corresponding to the circular symmetric sys- 
tem. In this case the density wave of a spherical wave is 
possible to appear as a stable ground state. 

In occurrence of the density wave instability there are 
two key ingredients, that is, the linear momentum depen- 
dence of V2d{Q) and the singular behavior of x {Q: ^/2) at 
q — 2ao^Fo- Thus the present mechanism of the den- 
sity wave instability in 2D dipolar Fermi gases is arisen 
from combined effects of the long-range nature of the 
dipole-dipole interaction and the Fermi surface effect in 
the 2D system. For 3D dipolar Fermi gases, the dipole- 
dipole interaction does not depend on the magnitude of 
momentum Q and the Fermi surface effect is less signif- 
icant. Thus the density wave phase may not appear in 
the 3D system. 

Figure [7] shows a phase diagram at zero temperature 
in terms of g and Oq. While the density wave instability 
is identified by Eq. (pS)) . the collapse instability is identi- 
fied with a negative value of the inverse compressibility. 



In this section, we apply our analysis for RPA stability 
into dipolar Fermi gases at finite temperatures and ob- 
tain the critical temperature of the density wave phase 
transition. 

To do so, we judge the stablity condition ([T71) with the 
density-density response function ([T8)) using the single 
particle energy calculated by the selfconsistent HF equa- 
tion (jlOp at a finite T without the variational ansatz. 
A critical temperature of the phase transition, Tc, is a 
highest temperature at which the condition 

^2<i(9,|)x(g,f;Te) = -i (27) 

is fulfilled. 

Figure [H] shows Tc in units of ideal gas Fermi temper- 
ature Tp = h^k^pi-j2mkB, as a function g for = 
(cross), 9a — 0.27r (asterisk), and 9^ = arccos(l/-\/3) 
(square). In Fig. [HI we also plot results by replacing e(fc) 
in Eq. ^ with the VM single particle energy of Eq. 
where A and ao are results of the corresponding system 
at zero temperature. The critical temperature with the 
HF single particle energy is higher than that with the 
VM single particle energy for same g and ^o- This in- 
dicates that the response function x(g,7r/2) at around 
q = kp(n /2) in HF approximation is larger than that in 
VM as shown in Fig. [51 This is a consequence of overesti- 
mate at fc > fci?(7r/2) for the VM single particle energy, as 
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FIG. 8: Critical temperature Tc with the HF single particle 
energy as a function of g for 6^0 = (cross), = 0.27r (aster- 
isk), and 6o — arccos(l/\/3) (square). The solid line, dashed 
line, dotted line are Tc with the VM single particle energy for 
^0 ~ 0, 0.27r, arccos(l/\/3), respectively. 
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FIG. 9: Density-density response function in units of m/h? 
for g = 0.928, So = a.rccos{l / y/i) , and T = 0.02rg in HF 
approximation (solid line) and VM (dashed line) . 



shown in Fig. [5l which leads to decrease of the response 
function via the denominator of integrand in Eq. (jl8p . 

It is well known that the charge density wave instabil- 
ity in ID conductor is due to the nesting property that 
is the increase of the number of lower energy excitations 
of particle-hole pairs near the Fermi surface [l^l- The 
fact that the curvature of the single particle energy at 
(/) = 7r/2 in HF approximation is lower than that in VM 
brings a large enhancement of the nesting property for a 
highly deformed Fermi surface. As shown in Fig. [Sj the 
difference between two results is remarkable for larger 
6*0- Thus we conclude that the deformation of the Fermi 
surface furthers the density wave phase transition. 

Finally, let us consider how high the critical temper- 
ature will be for a realistic system. For ''"K-^^Rb polar 
molecules with the electric dipole moment of 0.566 De- 
bye observed by JILA group [3], to = 127 a.m.u, and the 
number density of 2.0 x 10^ cm~^, we have g = 0.409. 
Figurelinireveals that Tc increases with decreasing 6*0 and 
can be a significant fraction of Tp . Such a high critical 
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FIG. 10: Critical temperature Tc as a function of 9o for g = 
0.409. 



temperature suggests that the density wave of fermionic 
polar molecules can be observed in future experiments. 



V. SUMMARY 

In the present paper we have studied a realization of 
the density wave phase transition in a 2D dipolar Fermi 
gas where the dipole moments of fermions are aligned 
by an orientable external field. Owing to judge the sta- 
bility of a normal gas, we have investigated RPA stabil- 
ity against density fluctuations of finite momentum. We 
showed that the density wave instability takes place in a 
broad region where the system is stable against collapse. 
We found that the critical temperature Tc of the phase 
transition can be achieved at a significant fraction of Tp 
for a realistic system of polar molecules. 

As discussed in Ref. pj|, the superfluid phase tran- 
sition of p-wave Cooper pairs is possible to realize for 
00 ^ arcsin(2/3) = 0.237r at T = 0. The combination of 
this and our conclusions offers a challenging quest for a 
coexisting phase of the superfluid and density wave or- 
ders in a 2D dipolar Fermi gas. This will be discussed in 
a future publication. 
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Appendix A: calculation of variational energy 

In this appendix, we obtain the variational energy at 
zero temperature under the assumption of the variational 
distribution function The total energy in HF ap- 

proximation is composed of the kinetic energy E^in and 
interaction energy Eint defined by Eqs. ^ and ([7]), re- 
spectively. 

First, the kinetic energy is easily obtained by 



_ n^nj^TT / I 2 

Ekin — — \ ^ + a 

Zm V 



Next, we calculate the interaction energy 
I r (fk r cPk' 



(Al) 



Ei, 



2 J (27r)2 J (27r)2 
^k%I{a;9a)C, 



f{k)f{k')V2d{k - k') 



where I{a;9o) is defined by Eq. (fT4l) and 

l-Tr/2 



c = 



noo poo /"TT/^ 

/ dxx^ / dyy / d(j)Q (1 ~ — y^ — 2xy cos ( 
Jo Jo Jo 



We define the function f{(f>) as 

f{(f>) = 1 - - 2/2 _ 2xy cos (f> > f.min 

where fmin = 1 - {x + y)^. If f,nin > then <d{f {(())) = 1 
for < < 7r/2, otherwise <d{f{(f))) = 1 for 0o < < vr/2 



and Q{f {(/))) = for < < 0o where the angle (j)o is 
determined by /(0o) = 0, that is 



00 — arccos 



With C = Ci - C2 where 



1 — a;2 _ y2 
2xy 



Ci = - dxx^ / dyy = 7r/30, 

^ Jo Jo 



and 



C2 = 



, /l-a:2-y2 

a;" / ay ?/ arccos I '■ — 

Ji-x \ 2a;?/ 



dx-^ 



= 7r/30 - 2/45, 
we obtain the interaction energy as 



Ei, 



32 h^n^ 



2d 



15 m 



gI{a;9o) 



(A2) 



Combining Eq. (|A1[) and Eq. (|A2[) . we obtain the varia- 
tional energy as 



rp 2a 

-C/fnt — 



vr / 1 



2 Va 



32 
15 



(A3) 
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